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ABSTRACT 

Single-epoch virial black hole (BH) mass estimators utilizing broad emission lines have been routinely 
applied to high-redshift quasars to estimate their BH masses. Depending on the redshift, different line 
estimators (Ha, H/3, MgllA2798, CivA1549) are often used with optical/near-infrared spectroscopy. Here 
we use a homogeneous sample of 60 intermediate-redshift (z ^ L5-2.2) SDSS quasars with optical and 
near-infrared spectra covering Civ through Ha to investigate the consistency between different single-epoch 
virial BH mass estimators. We critically compare restframe UV line estimators (CivA1549, Ciii]A1908 
and Mg II A2798) with optical estimators (H/? and Ha) in terms of correlations between line widths and 
between continuum/line luminosities, for the high-luminosity regime (L5100 > 10"*^"* ergs"') probed by our 
sample. The continuum luminosities of L1350 and L3000, and the broad line luminosities are well correlated 
with L5100, reflecting the homogeneity of quasar spectra in the restframe UV-optical, among which L1350 and 
the line luminosities for Civ and Cm] have the largest scatter in the correlation with L^ia). We found that 
the Mg II FWHM correlates well with the FWHMs of the Balmer lines, and that the Mg II line estimator 
can be calibrated to yield consistent virial mass estimates with those based on the H/3/Ha estimators, thus 
extending earlier results on less luminous objects. The Civ FWHM is poorly correlated with the Balmer 
line FWHMs, and the scatter between the Civ and H/3 FWHMs consists of an irreducible part (^ 0.12 dex), 
and a part that correlates with the blueshift of the Civ centroid relative to that of H/3, similar to earlier 
studies comparing Civ with Mg ll. The Cm] FWHM is found to correlate with the Civ FWHM, and hence 
is also poorly correlated with the H/3 FWHM. While the Civ and Cm] lines can be calibrated to yield 
consistent virial mass estimates as H/3 on average, the scatter is substantially larger than Mg II, and the usage 
of Civ/Ciii] FWHM in the mass estimators does not improve the agreement with the H/3 estimator. We 
discuss controversial claims in the literature on the correlation between Civ and H/3 FWHMs, and suggest 
that the reported correlation is either the result based on small samples or only valid for low-luminosity objects. 

Based in part on observations obtained with the 6.5 m Magellan-Baade telescope located at Las Campanas 
Observatory, Chile, and with the Apache Point Observatory 3.5 m telescope, which is owned and operated by 
the Astrophysical Research Consortium. 

Subject headings: black hole physics — galaxies: active — quasars: general 



1. INTRODUCTION 

Knowing the mass of active supermassive black holes 
(SMBHs) is of fundamental importance to understanding 
many physical processes associated with the black hole, as 
well as the assembly history of the SMBH population across 
cosmic time. Ove r the past several decades, reverberation 
mapp ing (RM. e.g..lBahcall et alJ [T97l iBlandford & McKed 
119821; lPetersonlll993h has proven to be a viable technique to 
measure (broad-line) AGN BH mass by providing an estimate 
of the broad line region (BLR) size R (e.g., IPeterson et al.l 
|2004|) . combined with the assumptions that the BLR dynam- 
ics is dominated by the central BH mass and that the widths 
of the broad emis sion lines V are related to the viri al velocity 
of the BLR (e.g.. IPibalfTgSOt IWandel et alJfT999h . The un- 
known geometry of the BLR is absorbed in a con stant virial 
coefficient /, which is calibrated ( e.g., lOnken e t al. 200'i; 
I Woo et al.ir2010tlGraham et al.ll201 Ih to bring the products of 
RV^/G into average agreement with those predicted from the 
local scaling relation between BH mass and bulge velocity 
dispersion (the M-a relation). 

An important result of RM studies is the discovery of a tight 
correlation between the BLR size and the continuum luminos- 
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ity of broad-hne AGNs (e.g., iKaspi etal J 120001: iBentz et alJ 
'2006'), i.e., the R-L relation, when plotted over a wide dy- 
namical range in AGN luminosity. This relation has led 
to the development of the so-called single-epoch virial BH 
mass estimators ("virial BH mass estimators" for short, e.g., 
Vestergaard 2002; McLure & Jarvis 2002; McLure & DunloEl 
I2004t IGreene & Holl2005t IVestergaard & Pete rson 2006), in 
which one measures the continuum (or line) luminosity and 
broad line width from single-epoch spectroscopy to derive 
a virial product as the BH mass estimate, with coefficients 
calibrated from a sample of ^ 40 local AGNs with RM 
masses (which are further tied to the predictions from the 
M-a relation). Various versions of single-epoch virial BH 
mass estimators have been developed since, based on dif- 
ferent broad lines and advocating different recipes for mea- 
suring luminosities and line widths (e.g ., McGill et al. 20081 
IWang et al."2009'; 'Vesterg aard & Osmeiif2009t iRafiee & Halll 
^lla; Shen et al. 2011). This empirical method, albeit 
rooted in the RM technique, is much less expensive than 
RM, and hence has been applied in numerous studies to es- 
timate quasar/ AGN BH masses, notably for large statistical 
samples (e.g., |W oo & Urrv"2002'; 'McLure & Dunlop"2004l; 
i Kollmeier et alTf2 006; Greene & Ho 2007; Vestergaard et afl 
l2008!; IShen et'aLi 2008. 2011). 

Despite the wide application of these virial BH mass esti- 
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mators, there are many statistical and systematic uncertain- 
ties of these estimates. First and foremost, all single-epoch 
mass estimators are bootstrapped from a sample of only ^ 40 
z < 0.4 RM AGNs (consisting of Seyfert 1 galaxies and 
several PG quasars), which is known to be unrepresentative 
of their high-lumino sity and high-redshift counterparts (e.g., 
[Richards et^l2011h . The statistics of RM AGNs need to be 
substantially improved to account for the diversity in BLR 
properties. Secondly, different versions of virial mass esti- 
mators have different systematics depending on the quality of 
the spectrum and the profile of the broad line, and there is 
currently no consensus as to which version is the best. Nev- 
ertheless, there are some general considerations on various 
estimators: 

• Which line to use: The commonly utilized pairs of line 
and luminosity in the restframe UV and optical are: Ha with 
Liia or L5100, H/3 with L5100, Mg II with L3000, and Civ with 
^^1350 (or L1450). Since the Balmer lines Ha and H/3 are the 
most studied lines in reverberation mapping and the R-L re- 
lation was originally measured for the Balmer line BLR radius 
and L5ioo(e.g., Kaspi et al. 2000; Bentz et al. 2006), it is rea- 
sonable to argue that the virial mass estimator based on the 
Balmer lines is the most reliable one. The width of the broad 
Ha is well correlated with that of the broad H/3 and there- 
fore it provides a good substitution in the absence of H/3 (e.g., 
I&eene & Ho 2005). 

The Mg II line has not been studied much in RM (cL. IWool 
I2OO8I) . and only in very f ew cases has a time-lag of Mg II been 
measured with RM (e.g. jReichert et al.lfl994tlMetzroth et all 
I2OO6I) : but the width of Mg II is shown to correlate with that 
of H/3 in single - epoch spectra (e.g. , Salviander et al. 2007; 
IShen et alJlIOOSt iMcGill etaUfMm. IWang et alJlIOOgh . sug- 
gesting that Mg II may be used as a substitution for H/3 in 
estimating virial BH masses. 

The Civ line is known to vary and time-lags have been mea- 
sured for Civ in several objects (e.g., Peterson et al. 2004; 
iKaspi et aP l2007h . although the sample is too small to de- 
rive a reliable R-L relation for Civ. However, the high- 
ionization Civ line differs from low-ionization lines such 
as Mg II and the Balmer lines in many ways (for a re- 
view, see Sulentic et al. 2000), most notably it shows a 
prominent blueshift with respect to the low-ionization lines 
(e.g., Gaskelli |1982|) . In addition, the Civ line is gen- 
erally more asymmetric than Mg II and the Balmer lines, 
and the width of Civ is poorly correlate d with those of 
Mgii and H/3 (e g.,|Bask in &Laor 2005; N etzer et al.|[2007l: 
IShen et al] l2008 f). The different properties of Civ suggest 
that Civ is probably more affected by a non-virial compo- 
nent such as arising from a radiatively- driven disk wind (e.g., 
iMurrav et all Il995i IProga et all 120001) . and would therefore 
be a biased virial mass estimator (e.g., JB askin & Laor 2005; 
Sulentic et al. 2007; Ne tzer et all I2007L IShen et al.1 I2008t 
Marziani & Sulentic..201 ll and references therein). However, 



since both the Balmer lines and Mg II move out of the optical 
bandpass at z > 2, it would be useful to improve the Civ esti- 
mator in order to measure BH masses at high redshift without 
the ne ed for n ear-IR spectroscopy. 

Shen et all ([2008) used a large sample of ~ 5000 SDSS 
quasars to show that the difference between the Civ and Mg II 
virial masses is correlated with the Civ-Mg II blueshift. On 
the other hand, lAssef et aTl (1201 Ih used a sample of ^ 10 
quasars with optical spectra covering Civ and near-IR spec- 
tra covering H/3/Ha to show that the difference between Civ 
and Balmer line virial masses is largely driven by their rest- 



frame UV-to-optical continuum luminosity ratio L1350/L5100, 
suggesting that much of the dispersion in their virial mass dif- 
ference is caused by the poor correlation between L5100 and 
L1350 rather than between their line widths. A larger sample is 
needed to test this result. 

• Line dispersion vs FWHM: The two common choices of 
line width are FWHM, and the second moment of the line 
(line dispersion, aum)- Both FWHM and crii„e have advan- 
tages and disadvantages. FWHM is easier to measure, less 
susceptible to noise in the wings and line blending than crune, 
but is more sensitive to the treatment of the narrow line re- 
moval. Arguably aijne is a b etter surrogate for the virial veloc- 
ity (e.g.. lCollin et aljr2006l) , although the evidence is not very 
strong. Since currently all the RM BH masses are computed 
using (Tiine.rms measured from the rms spectra ("Pe terson et al.l 
120041) . ideally one would like to use an^e, albeit not mea- 
sured from rms spectra, in single-epoch virial mass esti- 
mators. In practice, however, aune measured from single- 
epoch spectra depends on the quality of the spectra, line pro- 
file, and specific treatment of deblending, and could differ 
signifi cantly from one observa t ion/analysis to another (e.g. , 
Denn eveFd] l2009t iFine et all IMot iRafiee & Halj l20Tlbl: 
Assef et al.ll2011h . Therefore in terms of readiness and re- 
peatabiUty, aum is less favorable than FWHM in single-epoch 
virial mass estimators. For these reasons, we will not utilize 
line dispersion in the current study. 

In this paper we investigate the reliability of the UV virial 
mass estimators (in particular Civ) compared with the H/? 
(or Ha) estimator with a carefully selected sample of quasars 
with good Civ to Mg 11 coverage in optical SDSS spectra and 
our own near-IR spectra covering H/3 and Ha. Our sample 
probes the high-luminosity regime (Lsioo > 10"*^"* ergs"', or 
i'boi ^ 2.5 X 10"*'' ergs"') of quasars, and thus such a study 
will provide confidence on estimating virial BH masses for the 
most luminous quasars (such as z > 6 quasars). Our sample 
is substantially larger than earlier samples in similar studies, 
which enables us to draw more statistically significant conclu- 
sions. We are interested in examining empirical correlations 
between line widths and continuum luminosities of two differ- 
ent lines, and any dependence of their virial mass difference 
on specific quasar properties. We describe our sample and 
follow-up near-IR observations in ^ The procedure of mea- 
suring spectral properties is detailed in §|3]and the results are 
presented in ^ We discuss the results in ^and conclude in 
^ Throughout this paper we adopt a flat ACDM cosmology 
with flA = 0.7, ilo = 0.3 and Ho = 70kms"i Mpc"' . 

2. DATA 

2.1. Sample Selection 

We sele ct our targets from the SPSS DR7 quasar catalog 
(e.g.. Schneider et al.1 l20T0t IShen et al.1 IMl for follow-up 
near-IR spectroscopy with the following two criteria: 

• redshift between 1.5 and 2.2 and avoiding redshift 
ranges where the H/3 and Ha lines fall in the telluric 
absorption bands in the near-infrared; 

• with good (S/N > 10) SDSS spectra covering Civ 
through Mg II, and no broad absorption features or un- 
usual continuum shapes. 

These criteria by design selects luminous quasars (bolo- 
metric luminosity Lboi > a few x 10"*^ ergs"') as our targets, 
but the resulting sample still covers a range of spectral di- 
versities such as the line width and velocity shift of each 
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broad lines. In addition, host contamination is generally neg- 
ligible for these objects, which greatly simplifies our model 
fits and interpretations. Our targets have a similar color ex- 
cess A{g - i) distribution (where A{g - i) is the deviation 
of g - i color from the median g - i color at each redshift; 
see [Richards et al.ll2003h as the underlying SDSS quasars in 
this redshift range, but do not include any dust-reddened ob- 
jects (defined as A(g-i) > 0.3). - 10% (5/49) of the tar- 
gets are radio-loud (R = /,y,6cm//i, 2500A ^ ^0, see lShen et alj 
l2011h based on the Faint Images qf the Radio Sky at Twenty- 
Centimeters (FIRST. IWhite et al.|[l997[) catalog, and the rest 
1 1 targets are not in the FIRST footprint as of July 16, 2008. 

2.2. Near-IR Spectroscopy 

We observed our targets during 2009-201 1 with TripleSpec 
dWilson et al. 20 04) on the ARC 3.5 m tel escope, and with the 
Folded-port InfraRed Echellette (FIRE, ISimcoe et all 120101) 
on the 6.5 m Magellan-Baade telescope. Table [T| summarizes 
our sample and follow-up observations. Below we describe 
the observations and data reduction for TripleSpec and FIRE 
data, respectively. 

2.2.1. ARC 3.5m/TripleSpec 

TripleSpec (IWilson et al.l 120041) is a near-IR spectrograph 
with simultaneous 0.95- 2.46 /im overage. We observed our 
targets during 2009-2011 semesters. The total exposure time 
varied from object to object due to different target brightness 
and observing conditions, but is typically 1-1.5 hr We used 
slits with widths of both 1.1" and 1.5" during the course of 
the observations, and the resulting spectral resolution is /? ^ 
2500 - 3500. The slit was positioned at the parallactic angle 
in the middle of the observation, and we performed standard 
ABBA dither patterns to aid sky subtraction. For each object 
we observed a nearby AOV star as flux and telluric standard 
immediately before or after observing the science target. 

We reduced the Triplespec data using the IDL-based 
pipeline APOTripleSpecTool, which is a modified ver- 
sion of_flie_S2e5ctool_package developed by Michael Gush- 
ing dCushing et alj|2004h . The reduction procedures include 
non-linearity correction, flat-fielding, wavelength calibration 
using OH sky lines (calibrated to vacuum wavelength), sky 
subtraction using adjacent exposures at nodding slit posi- 
tions, cosmic -ray rejection, optimal extraction of 1-D spectra 
(iHornd 19861) . combining individual exposures, merging mul- 
tiple echelle orders, and heliocentric corrections. We used the 
AOV-star observations for relative flux calibration a nd tell uric 
correction following the technique of lVaccaetal.l (l2003l) us- 
ing the xtellcor routine contained in the Spextool package 
(Gushing et al. 2004). We tied the absolute flux calibration 
to the Two Micron All Sky Survey (2MASS) (Skrutsk ie etal] 
I2OO6) //-band magnitudes using synthetic magnitudes com- 
puted from our spe ctrum with the 2M ASS relative spectral 
response curves in iGohen et al.l (l2003l) . This absolute flux 
calibration neglects the continuum variability between the 
2MASS and (spectroscopic) SDSS epochs, which is typi- 
cally at the level of ^ 0.1 mag for average SDSS quasars 
(e.g.. ISesar etal .1120071: iMacLeod et al.1 1201 ll) . It also ne- 
glects possible line shape variability of quasars between the 
two epochs of SDSS and near-IR observations, but this vari- 
ation is likely negligible (ctfwhm < 0.05 dex) base d on re- 
peated spectrosc opy of the same objects (e.g., Wilhi te et alJ 
l20()7HpSce^a^l2012^ . 

Finally, for 14 targets we have a second observation on a 
different night. We combined these repeated observations us- 



ing the inverse-variance weighted mean of the two observa- 
tions. 

2.2.2. Magellan/FIRE 

FIRE JSimcoe et al.ll2010l) is a near-IR echelle spectrome- 
ter covering the full 0.8-2.5 /im band. We observed 20 targets 
during the nights of April 25-26, 2011, and another two tar- 
gets on the nights of July 12-13, 2011. We used the 0.6" 
slit width in Echelle mode, which offers a spectral resolution 
of R ^ 6000 (50kms"'). Typical total exposure times were 
45 min per target but varied from object to object. We ob- 
served our targets at the parallactic angle, and for each target 
we observed a nearby AOV star for flux and telluric standard. 

We reduced the FIRE data using the IDL-based pipeline 
"FIREHOSE" developed by Robert Simcoe et al ^. The re- 
duction procedures are similar to those for the TripleSpec data 
with the exception of sky subtraction. Instead of subtracting 
adjacent nodding exposures, sky subtraction was performed 
using a B-spline model of the sky dir ectly constructe d for each 
exposure following the technique of lKelsonl (l2003l) . 

The FIRE spectra have a substantial spectral overlap with 
the SDSS spectra. Therefore we used the common part with 
the SDSS spectrum to normalize the FIRE spectral flux den- 
sity. As for our TripleSpec data, we neglect variations in line 
shape between the SDSS and FIRE spectroscopic epoches. 

3. SPECTRAL MEASUREMENTS 

To derive line width and continuum luminosities used in 
single-epoch virial mass estimators, we perform spectral fits 
to the optical and ne ar-IR spectra, as commonly adopted 
in the literature (e.g., Greene & Ho 2005; Salviande ret al.l 
l2007t IShen et al. 2008; Wang et al. 2009). Spectral fits with 
some functional form have certain advantage of being less 
susceptible to noise than direct spectral measurements, al- 
though sometimes there are still ambiguities in decomposing 
the spectrum into different components. 

Here we perform least- x^ global fits to the combined opti- 
cal and near-IR spectra for the same object. Such global fits 
were not possi ble for object s with limited wavelength cov- 
erage (e.g., She n et al]|2011l) . Each combined spectrum was 
de-reddened for Galactic extinction using the iGardelh et aP 
( 19 89) Milky Way redde ning law and E{B - V) derived from 
the ISchlegel et alj (Il998h dust map. The spectrum was then 
shifted to restframe us ing the improved redshifts provided by 
iHewett & Midi (12010 ) for SDSS quasars, where the spectral 
fit was performed. For each object we masked out narrow ab- 
sorption line features imprinted on the spectrum, which will 
bias the continuum and emission Une fits. 

3.1. Pseudo-Continuum Fit 

We first fit a pseudo-continuum model to account for 
the power-law (PL) continuum, Fe II emission and Balmer 
continuum underneath the broad emission lines of inter- 
est. All components were fit simultaneously. Templates 
for Fe II and Fe III emission have been constructed from 
the s pectrum of the narrow-line Seyfert 1 galaxy, I Zw 1 
(e.g., 'Bor oson & Greeiil 119921; IVestergaard & Wilkesl l2001l 
Tsuzuki et al] |2006l) . In this work we do not include addi- 
tional Fe III emission in the fits as we found this component 
is poorly constrained ( e.g., Greene et al. 2010 ) . For t he UV Fe 
II template, we use the lVestergaard & Wilkesl (12001 1) template 

' http://web.mit.edu/~rsimcoe/www/FIRE/ob_data.htm 
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Sample Summary 












Object Name 


RA (J2000) 


DEC (J2000) 


Plate 


Fiber 


MJD 


ZHW 


ipsF 


-'2MASS 


W2MASS 


^s.2MASS 


NIR Obs. 


Obs. UT 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


J0029-0956 


00 29 48.04 


-09 56 39.4 


0653 


640 


52145 


1.618 


17.672 


16.728 


15.747 


15.622 


TSPEC 


100102/101128 


J004 1-0947 


00 41 49.64 


-09 47 05.0 


0655 


172 


52162 


1.629 


16.966 


16.201 


15.680 


15.535 


TSPEC 


100102/101128 


J0147+1332 


01 47 05.42 


+13 32 10.0 


0429 


145 


51820 


1.595 


17.115 


16.194 


15.473 


15.474 


TSPEC 


090909/091107 


J0149+1501 


01 49 44.43 


+15 01 06.6 


0429 


575 


51820 


2.073 


17.275 


16.565 


15.998 


15.243 


TSPEC 


090909/101128 


JO 157-0048 


01 57 33.87 


-00 48 24.4 


0403 


213 


51871 


1.551 


18.164 


16.797 


16.553 


0.000 


TSPEC 


091107/101128 


J0200+1223 


02 00 44.50 


+12 23 19.1 


0427 


219 


51900 


1.654 


17.811 


16.573 


16.125 


0.000 


TSPEC 


100102/101128 


J0358-0540 


03 58 56.73 


-05 40 23.4 


0464 


499 


51908 


1.506 


18.258 


17.572 


16.297 


0.000 


TSPEC 


100102/101128 


J04 12-06 12 


04 12 55.16 


-06 12 10.3 


0465 


037 


51910 


1.691 


17.322 


16.306 


16.077 


15.353 


TSPEC 


100102/101128 


J0740+2814 


07 40 29.82 


+28 14 58.5 


0888 


545 


52339 


1.545 


17.445 


16.426 


15.689 


15.482 


TSPEC 


091108 


J08 12+0757 


08 12 27.19 


+07 57 32.9 


2570 


026 


54081 


1.574 


17.404 


16.658 


15.995 


16.031 


TSPEC 


101202 


J08 13+2545 


08 13 31.28 


+25 45 03.0 


1266 


219 


52709 


1.513 


15.385 


14.085 


13.271 


13.056 


TSPEC 


091108 


J0813+1522 


08 13 44.15 


+15 22 21.5 


2270 


439 


53714 


1.545 


17.541 


16.472 


15.805 


0.000 


TSPEC 


101122 


J0821+5712 


08 2146.22 


+57 12 26.0 


1872 


615 


53386 


1.546 


16.868 


15.943 


15.027 


15.031 


TSPEC 


091108/100104 


J0838+2611 


08 38 50.15 


+26 1105.4 


1930 


492 


53347 


1.618 


16.098 


15.211 


14.424 


14.288 


TSPEC 


091108 


J0844+2826 


08 44 51.91 


+28 26 07.5 


1588 


179 


52965 


1.574 


18.006 


17.026 


16.147 


15.798 


TSPEC 


101202 


J0855+0029 


08 55 43.26 


+00 29 08.5 


0468 


HI 


51912 


1.525 


17.952 


16.829 


16.545 


15.668 


FIRE 


110426 


J09 17+0436 


09 17 54.44 


+04 36 52.1 


0991 


284 


52707 


1.587 


18.543 


0.000 


0.000 


0.000 


FIRE 


110427 


J0933+1413 


09 33 18.49 


+14 13 40.1 


2580 


347 


54092 


1.561 


17.465 


16.520 


15.540 


15.448 


TSPEC 


100126 


J0941+0443 


09 41 26.49 


+04 43 28.7 


0570 


379 


52266 


1.567 


17.824 


16.954 


16.084 


15.706 


FIRE 


1 10427 


J0949+1751 


09 49 13.05 


+17 5155.9 


2370 


184 


53764 


1.675 


17.143 


16.137 


15.626 


15.332 


TSPEC 


100126 


J 1004+4231 


10 04 01.27 


+42 3123.1 


1217 


573 


52672 


1.666 


16.764 


15.795 


15.376 


15.080 


TSPEC 


100104 


J 1009+0230 


10 09 30.51 


+02 30 52.4 


0502 


429 


51957 


1.557 


18.556 


17.310 


0.000 


0.000 


FIRE 


110426 


J1014+5213 


10 14 47.54 


+52 13 20.2 


0904 


259 


52381 


1.552 


17.334 


16.705 


15.929 


15.752 


TSPEC 


110124 


J1015+1230 


10 15 04.75 


+12 30 22.2 


1745 


148 


53061 


1.703 


17.400 


16.374 


15.989 


15.658 


TSPEC 


110124 


J1046+1128 


10 46 03.22 


+11 28 28.1 


1601 


193 


53115 


1.607 


17.784 


0.000 


0.000 


0.000 


FIRE 


110426 


J1049+1432 


10 49 10.31 


+14 32 27.1 


1749 


571 


53357 


1.540 


17.740 


16.813 


15.590 


15.458 


TSPEC 


100126 


J 1059+0909 


10 59 51.05 


+09 09 05.7 


1220 


231 


52723 


1.690 


16.771 


15.620 


15.094 


14.411 


TSPEC 


110124 


Jl 102+3947 


11 02 40.16 


+39 47 30.1 


1437 


205 


53046 


1.664 


17.605 


16.563 


16.153 


15.811 


TSPEC 


1 10222 


Jl 119+2332 


11 19 49.30 


+23 32 49.1 


2493 


077 


54115 


1.626 


17.338 


16.230 


15.465 


15.322 


TSPEC 


110124 


Jl 125+0001 


11 25 42.29 


+00 0101.3 


0280 


077 
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Note. — Summary of the sample of SDSS 


quasars 


for which we have conducted near-infrared spectroscopy. 
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Fig. 1 . — An example of our model fits to the combined optical and near-IR spectra (J0941+0443). The top panel shows the global fit of the pseudo-continuum, 
where the brown line is the power-law continuum, the blue line is the Fe II template fit, the cyan line is the Balmer continuum model, and the red line is the 
combined pseudo-continuum model to be subtracted off. The line segments near the top indicate the wavelength windows used for the pseudo-continuum fit. The 
bottom panels show the emission line fits to CIV through Ha, where the cyan fines are the model narrow fine emission, the green lines are the model broad line 
emission, and the red fines are the combined model line profiles. For CIII] we also show the modeled AlIII and Silll] emission in magenta. 
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(1000-3090A). ISalvianderefai] (l2007b modified this tem- 
plate by extrapolating below the Mg II line, and we use their 
template for the 2200-3090 A region; we augment the 3090- 
3500 A region using the template derived by iTsuzuki et alj 
( 120061) . For the opti cal Fe II template (3686-7 484 A) we use 
the one provided bv lBoroson & GreenI (11992 1). The PL con- 
tinuum model has two free parameters, the normalization and 
the PL slope. The UV and optical Fe II templates are fitted 
independently, each has three free parameters, the normaliza- 
tion factor, the velocity dispersion (to be convolved with the 
template), and the wavelength shift of the template. The Fe II 
templates are only used as an approximation to remove signif- 
icant iron emission, and we found that they did a reasonably 
good job. However, we are not concerned with the properties 
of the iron emission in this work, and thus we do not interpret 
the physical meanings of the velocity dispersion and wave- 
length shift of the Fe II templates. 

For the B aimer continuum we follow the empirical model 
by iGrandil (1982) as composed of partially optically thick 
clouds with an effect i ve temperature (e .g., iDietrich et al.l 
I200I IWangeFaDIIOOgt lOreene et al.ll2010h : 



/Bc(A)=ABA(A,r,)(l-e-^-); A < Abe 



(1) 



where Te is the effective temperature, Abe = 3646 A is the 
Balmer edge, t\ = tbe(A/Abe)"^ is the optical depth with tbe 
the optical depth at Abe, A is the normalization factor, and 
Ba(A, Te) is the Planck function at temperature Tf.. During the 
continuum fits, we have three free parameters, 1 x 10^ < Tg < 
5 X lO-'K, 0.1 < Tbe < 2, and A > 0. 

Note that for Umited wavelength fitting range (i.e., A < 
3646 A), the Balmer continuum cannot be well constrained 
and is degenerate with t he power-law and Fe II components 
(e.g.. IWang et al.l l2009l) . and is generally not fitted (e.g., 
IShen etal]l2011h . In the case of global fits, a single power- 
law continuum is required to simultaneously fit the region 
from Civ to Ha, providing some additional constraints on 
the Balmer continuum; however even in this case, the Balmer 
continuum may still be poorly constrained in a few cases, 
which will lead to uncertainties in the power-law continuum 
luminosity estimates. Nevertheless, the isolation of broad 
emission lines is not affected much by including or exclud- 
ing the Balmer continuum model. 

We fit the pseudo-continuum model to a set of continuum 
windows free of strong emission lines (except for Fe II): 1 350- 
1360A, 1445-1465A, 1700-1705A, 2155-2400A, 2480- 
2675 A, 2925-3500 A, 4200-4230 A, 4435-4700 A, 5100- 
5535 A, 6000-6250 A, 6800-7000 A. We try fitting both with 
and without the Balmer continuum component and adopt the 
fit with the lower reduced ^ value; usually adding the Balmer 
continuum improves the global fit. 

In a few cases (^ 5 objects) we found that this global 
pseudo-continuum model does not fit the Civ-Ciii] region 
well, which is likely caused by intrinsic reddening in these 
systems, ill-determined Balmer continuum strength, or mis- 
matched iron template. For these objects we perform lo- 
cal (A < 2165 A with the same continuum windows defined 
above) continuum fits around the Civ and Cm] regions with- 
out the Balmer continuum and Fe II emission (i.e., only with 
the PL component), in order to get better measurements for 
Civ and Cm]. The PL continuum model is then used to esti- 
mate the monochromatic continuum luminosity L\ = Xf\ at 
5 100 A, 3000 A and 1350 A. Although some of our targets 
do not have spectral coverage of the restframe 1350 A due 



to their relatively lower redshift, the global PL component is 
well constrained so the extrapolation of the model continuum 
to 1350 A is not a problem. 

3.2. Emission Line Fits 

Once we have constructed the pseudo-continuum model, 
we subtract it from the original spectrum, leaving the 
emission-line spectrum. We then fit the Ha, H/3, Mg II, Cm], 
Civ broad line complexes simultaneously with mixtures of 
Gaussians (in logarithmic wavelength), as detailed below: 

• Ha: we fit the wavelength range 6400-6800 A. We use 
up to 3 Gaussians for the broad Ha component, 1 Gaus- 
sian for the narrow Ha component, 2 Gaussians for the 
[N II] A6548 and [N 11] A6584 naiTow fines, and 2 Gaus- 
sians for the [S 11] A6717 and [S 11] A6731 narrow lines. 
Since the narrow [N 11] lines are underneath the broad 
Ha profile, we tie their flux ratio to be /6584//6548 = 3 
to reduce ambiguities in decomposing the Ha complex. 

• H/3: we fit the wavelength range 4700-5 100 A. We 
use up to 3 Gaussians for the broad H/3 component 
and 1 Gaussian for the narrow H/3 component. We 
use 2 Gaussians for the [Oiii] A4959 and [Oiii] A5007 
narrow lines. Given the quality of the near-IR spec- 
tra, we decided to only fit single Gaussians to the 
[O III] AA4959,5007 lines, and we tie the flux ratio of 
the [O III] doublet to be /5007//4959 = 3. 

• Mgll: we fit the wavelength range 2700-2900 A. We 
use up to 3 Gaussians for the broad Mg II component 
and 1 Gaussian for the narrow Mg II component. We do 
not try to fit the Mg II lines as a doublet, as the line split- 
ting is small enough not to affect the broad line width 
measurements, and the spectral quality is usually inade- 
quate for fitting a doublet to the narrow Mg II emission. 

• Cm]: we fit the wavelength range 1820- 1970 A. We 
use up to 2 Gaussians for the broad Cm] component 
and 1 Gaussian for the narrow Cm] component. We 
use two additional Gaussians for the Sim] A 1892 and 
Aim A 1857 lines adjacent to Cm]. To reduce ambigui- 
ties in decomposing the Cm] complex, we tie the cen- 
troids of the two Gaussians for the broad Cm] i.e., the 
broad Cm] profile is forced to be symmetric; we also 
tie the velocity offsets of Sim] and Aim to their relative 
laboratory velocity offset. 

• Civ: we fit the wavelength range 1500-1600 A. We use 
up to 3 Gaussians for the broad Civ component and 1 
Gaussian for the narrow Civ component. We do not fit 
the 1640 A Hen feature as its contribution blueward of 
1600 A is negligible and will not bias the Civ line fit. 

• During the line fitting, all narrow line components are 
constrained to have the same velocity offset and line 
width. We also impose an upper limit of 1200 km s"' 
for the FWHM of the narrow line component"*. 

'' This upper limit is slightly larger than the values used in some studies 
(typically ~ 750- 1000 kms"' ). For luminous SDSS quasars, [O III] FWHM 
values exceeding ~ lOOOkms"' are often seen (e.g.. IShen et aUlSoTl^ . We 
hereby adopt the 1200 kms"^' upper limit for the narrow line width. 



VIRIAL BH MASS ESTIMATORS 



While the presence of narrow Hne components for the 
Balmer lines is beyond doubt, the relative contribution from 
narrow line components for the UV lines is less certain. For 
Mg II, there is clear evidence that a narrow line component 
is present at lea st in some quasars (e.g., |S hen et al. 20081 
1201 U l Wangetal. 200 9). For Civ, the lVestergaard & PetersonI 
( |2006|) virial mass calibration uses the FWHM from the whole 
line profile, while some argue that a narrow line component 
shoul d be subtracted for Civ as well (e.g., iBaskin & Laoij 
12005^). The presence of narrow emission lines in the restframe 
optical spectra is essential to provide constraints on the nar- 
row line contribution for the UV lines. We will measure the 
Civ line width both with and without narrow line subtraction 
and test if it is necessary to remove narrow line emission for 
Civ. 

3.3. Measurement Uncertainties 

It is important to quantify the uncertainties in our spec- 
tral measurements. The nature of the non-linear model and 
multi-component fits introduces ambiguities in decomposi- 
tion, and the resulting uncertainties are usually larger than 
those estimated from the parameter co-variance matrix of the 
least- x^ fits. We estimate the uncertainties in mea sured spec- 
tral qu antities using a Monte-Carlo approach as in lShen et alJ 
(1201 lb . For each object we generate 50 random realizations 
of mock spectra by adding Gaussian noise to the original 
spectrum at each pixel using the spectral error array. Tech- 
nically speaking, these mock spectra are not the exact alterna- 
tive realizations of the original spectrum since the errors were 
added twice, but they are only slightly degraded realizations 
and provide a good approximation to capture the wavelength- 
dependent noise properties. We fit each mock spectrum with 
the same fitting procedure described above and derive the dis- 
tribution of each measured spectral quantity (such as FWHM, 
velocity offset, etc). We then take the semi-quartile of the 
68% range of the distribution as the nominal uncertainty of 
the measured quantity. This approach takes into account the 
statistical uncertainties due to flux errors, and systematic un- 
certainties due to ambiguities in decomposing multiple com- 
ponents. 

Fig. [T] shows an example of our global fits, and we tabu- 
late the measured quantities in Table |2] Although with dif- 
ferent fitting recipes, the continuum and emission line mea- 
surements are c onsistent with t he measurements with SDSS 
spectra alone (Shen et ani2011l) . and the largest discrepancy 
occurs for LMgii.broad and L3000: logLsooo is systematically 
smaller by ^ 0.12 dex, and logLMgii.bioad is sys tematically 
larger by ^ 0.067 dex, than the measurements in IShen et al.l 
(1201 Ih . This is largely caused by the additional Balmer con- 
tinuum model in the spectral fits. For simplicity, from now on 
we will by default refer to the broad line component when we 
mention the FWHM or luminosity of a particular line unless 
stated otherwise. 

4. RESULTS 

We now proceed to examine correlations between contin- 
uum (line) luminosities and between line widths for different 
lines, as well as their virial products. 

4.1. Luminosity Correlations 
In Fig.|2]we compare different luminosities with the con- 
tinuum luminosity at 5 100 A, and we list the slopes from 
the bisector linear regression fits using the BCES estima- 
tor dAkritas & BershadvllT996h in Table |3] Our objects all 
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Fig. 2. — Correlations of different luminosity indicators with L5100 for 
our sample. Each correlation has been shifted vertically for clarity without 
changing the scatter in the coiTelation. Upper: Correlations between L5100 
and the three most-frequently used altemative luminosity indicators for the 
BLR size, L^ooo- ^1350 ^nd Lna. broad- Bottom: Correlations between L5100 
and broad Une luminosities. In both panels the solid lines are the bisector 
linear regression results using the BCES estimator (e.g., Akiitas & B ershadvl 
I199fl) . and the dashed lines indicate a linear correlation of unity slope. 



have luminosity L5 100 > 10"*^"* ergs"', and therefo re contami- 
nation from host starlight is mostly negligible (e.g. JShen et al.l 
l2011h . For UV estimators, L3000 and L1350 (or Li4go) are often 
used in replacement of Lsioo- In addition, the luminosity of 
the b road Ha line Lna is also used in pair with Ha line width 
(e.g.. lGreene & Holl2005h . Since the original 7J-L relation is 
calibrated against Lsioo, a good correlation with Lsioo is re- 
quired to produce a reasonable estimate of the BLR size with 
alternative luminosity indicators. As shown in the top panel 
of Fig. |2] all three luminosity indicators are correlated with 
L5100, where Lhu and L3000 are better correlated with L5100 
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than Li35o^. The slopes of these correlations are very close to 
unity, which means the ratios of these luminosities to L5100 al- 
most do not depends on luminosity over this luminosity range. 

In cases where the continuum is too faint to detect or con- 
taminated by host starlight or emission from a relativistic jet, 
an alt ernative route is to use the lurninosity of the broad line s 
(e.g.. IWu et aHIIOOl iGreene & H0II2OO5E IShen et aUlIoTTh . 
The bottom panel of Fig. |2] shows correlations of different 
line luminosities with Lsioo- Again, the H/3 and Mgll line 
luminosities seem to correlate with L5100 with lower scatter 
than Cm] and Civ. Interestingly, the scatter in the Lciv-^sioo 
relation is comparable to that in the Logo-Lsioo relation, sug- 
gesting that using the Civ line luminosity will not degrade the 
mass estimates much than using Logo. The best-fit slopes of 
these correlations are slightly different from unity, indicating 
a possible mild luminosity dependence of the ratios of these 
luminosities to Lgioo over this luminosity range. 

We note that these luminosity correlations are not predom- 
inately caused by the common distance of each object. In 
fact, the dynamical range resulting from luminosity distances 
is only 0.4 dex given the limited redshift range of our objects, 
while the entire luminosity span is 1 .5 dex. These luminosity 
correlations justify the usage of alternative luminosity indi- 
cators in various virial mass estimators. The scatter in the 
correlation between alternative luminosity indicator and L5100 
will be one source of the scatter in virial mass estimates when 
compared with those based on H/3 width and Lgioo- 

4.2. Line Width Correlations 

While it is still debated whether or not a narrow line com- 
ponent for Civ needs to be subtracted for high-redshift broad- 
line quasars (see d iscussions in, e.g., iBachev et al.l l2004i 
ISulentic et al.ll2007h . it is clear that narrow Civ emission doe s 
exist, as seen in some type 2 quasars (e.g.. lStern et al.ll2002l) . 
The [O III] coverage in our near-IR spectra makes it possible 
to constrain the strength of the narrow Civ emission by fix- 
ing its line width and velocity offset to those of the narrow 
[Olll] lines. We have measured the width of Civ with and 
without the subtraction of a possible narrow line component. 
In Fig. [3] we compare the resulting Civ FWHM with the two 
methods. The two objects with large error bars (J1009+0230 
and J1542+1112) have associated absorption, which causes 
some ambiguities in decomposing the Civ line in our Monte 
Carlo mock spectra and therefore leads to large uncertain- 
ties (see ®. We found that the narrow line contribution to 
Civ is generally weak for objects in our sample, and only 
in 2 objects (Jl 119+2332 and J1710+6023) the narrow Civ 
component is strong enough to make a large difference in the 
line width measurement. From now on we use the Civ line 
width with narrow line subtraction. We note, however, that 
our quasars are luminous, and the relative strength of the nar- 
row Civ emission may be larger fo r lower lu minosity objects 
(see, e.g. jBachev et alJl20 04: Sule ntic et al.l l2007). 

In the bottom panel of Fig. [3] we compare the Civ FWHM 
and Cm] FWHM. The measurement errors are typically larger 
for Cm] due to the ambiguity of decomposing the Cm] com- 
plex, but a correlation is still seen between the FWHMs of 
Civ and Cm]. This is intriguing because Cm] does not show 
as large a blueshift relative to the low-ionization lines as does 

^ The measured L3000 with the Balmer continuum component in the fit 
is on average smaller by ^ 0. 12 dex than that without fitting the Balmer 
continuum. However, both measurements of L3000 are tightly correlated with 
^5100- 
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Fig. 3. — Upper: Comparison between the two methods of measuiing the 
CIV FWHM, i.e., with and without subtracting a narrow line component. The 
two methods yield similar results for the majority of our objects, indicating 
that the naiTow line contribution is generally negligible for CIV for luminous 
quasars with L5100 > 10'*''* ergs"'. Bottom: Comparison between the CIV 
and cm] FWHMs. The uncertainties associated with the CIII] FWHM mea- 
surements are typically large due to ambiguities in decomposing the CIII] 
complex, but a general coiTelation is seen between the two FWHMs. 



Civ (e.g.. i Richards et al]|2011l Shen et al., in preparation) 
when the contributions from Sim] and Aim are removed. In 
Fig.lHwe plot the FWHM against the velocity offset relative to 
the broad H/3 line, for Civ and Cm] respectively. In both cases 
a significant positive correlation is detected, i.e., FWHM in- 
creases with blueshift. Such a trend is already known when 
compai-ing Civ and Mgn (e.g.. IShen et al.ll2008l [201 Ih . and 
now it is confirmed when comparing Civ directly with H/3. 
However, similar trends were not found for the FWHM of 
Ha, H/3 or Mg n against their velocity shift relative to [O III]. 
We also tested if there is any correlation between FWHM and 
other properties (continuum luminosity, color, line asymme- 
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TABLE 2 
Spectral Measurements 



FoiTnat Units 



Description 



objname 


AlO 


- 


logil350 


F6.3 


ergs"' 


ErrlogLi35o 


F6.3 


ergs"' 


logLsooo 


F6.3 


ergs"' 


En-IogL3ooo 


F6.3 


ergs"' 


log is 100 


F6.3 


ergs"' 


En- log L5 100 


F6.3 


ergs"' 


logLciv 


F6.3 


ergs"' 


ErrlogLciv 


F6.3 


ergs"' 


logi-cill] 


F6.3 


ergs"' 


ErrlogLciii] 


F6.3 


ergs"' 


logi-Mgll 


F6.3 


ergs"' 


EnlogLMgii 


F6.3 


ergs"' 


logiH/3 


F6.3 


ergs"' 


ErrlogLH/3 


F6.3 


ergs"' 


lOgiHa 


F6.3 


ergs"' 


EnlogLHQ 


F6.3 


ergs"' 


FWHMciv 


15 


kms" 


En- FWHMciv 


15 


kms"' 


FWHMciii] 


15 


kms"' 


EnFWHMcmj 


15 


kms" 


FWHMMgii 


15 


kms"' 


En FWHMMaii 


15 


kms"' 


FWHMhs 


15 


kms" 


En FWHMh;3 


15 


kms"' 


FWHMHa 


15 


kms"' 


EnFWHMHc, 


15 


kms" 


VciV-H;3 


15 


kms"' 


En Vciv-HS 


15 


kms"' 


VciV-AlIII 


15 


kms"' 


En Vcrv-Aiiii 


15 


kms"' 


Vcin]-H/3 


15 


kms"' 


En Vcm]-H^ 


15 


kms"' 


VMgII-[OIII] 


15 


kms" 


En yMgii-[om] 


15 


kms"' 


VH/3-[0in] 


15 


kms"' 


En yH;3-[omj 


15 


kms"' 


VHa-[Om] 


15 


kms" 


En yHa-[OIIIJ 


15 


kms"' 


CIV AS 


F4.2 


- 



Object Name 
continuum luminosity at restframe 1 350 A 

measurement enor in log L 1350 
continuum luminosity at restframe 3000 A 

measurement enor in log L3000 

continuum luminosity at restframe 5100 A 

measurement enor in log L5 100 

luminosity of the broad CIV line 

measurement enor in log Lciv 
luminosity of the broad CIII] line 

measurement enor in log Lcni] 

luminosity of the broad Mg II line 

measurement enor in log LMgii 

luminosity of the broad H/3 line 

measurement enor in log L^^ 

luminosity of the broad Ha line 

measurement enor in log Lua 

FWHM of the broad CIV line 

measurement enor in FWHMciv 

FWHM of the broad CIII] line 

measurement enor in FWHMciii] 

FWHM of the broad Mg II Une 
measurement enor in FWHMMgii 

FWHM of the broad H/3 line 
measurement enor in FWHMh^ 

FWHM of the broad Ha line 

measurement enor in FWHMhq 

blueshift of the broad CIV centroid w.rt. the broad H/3 centroid 

measurement enor in Vciv-H/3 
blueshift of the broad CIV centroid w.rt. the broad AlIII centroid 

measurement enor in Vciv-Allll 

blueshift of the broad CIII] centroid w.rt. the broad H/3 centroid 

measurement eiTor in Vcinj-H/3 

blueshift of the broad Mg II centroid w.rt. the nanow [O III] centroid 

measurement error in VMgii_[oiiij 
blueshift of the broad H/3 centroid w.rt. the naiTow [O III] centroid 

measurement enor in Vh/S-ioiii] 

blueshift of the broad Ho centroid w.rt. the nanow [O III] centroid 

measurement enor in VHa-[Olll] 

Asymmetry parameter of the broad CIV line 



Note. — Format of the tabulated spectral measurements. The full table is available in the online version. 



TABLE 3 
BCES Bisector Slopes for 
Luminosity Correlations 



VS L5100 


a 


(To, 


Scatter 


^-1350 


1.044 


0.099 


0.13 


^■3000 


0.979 


0.036 


0.05 


'^H a. broad 


I.OIO 


0.042 


0.07 


^H/S.broad 


1.251 


0.067 


O.ll 


^Mgll.broad 


0.861 


0.070 


O.ll 


icilll.broad 


0.924 


0.099 


0.16 


iciV.broad 


0.950 


0.085 


0.14 



Note. — a and Ca are the slope 
and uncertainty (Icr) in slope from the 
bisector linear regression fit of each lu- 
minosity against Lsioo- The last col- 
umn lists the scatter perpendicular to 
the best-fit linear relation (dominated 
by intrinsic scatter rather than measure- 
ment errors). 



try) for all five lines and found none of them is significant. 
In Fig. |5] we plot different line widths against the broad 



H/3 FWHM. We have suppressed measurement errors in these 
plots for clarity. In addition to the traditional FWHM, we also 
measure the full-width-at-third-maximum (FWTM) and full- 
width-at-quarter-maximum (FWQM) as alternative line width 
indicators. Consistent with earlier studies, we see strong 
correlations among t he widths of Ha, H/ 3 and Mg II (e.g. , 
Greene & H o 2005; Salviander et aUlIOOTI: IShen et al.ll200l 
Wang et al.|[2009i) . On the other hand, both Cm] and Civ line 
widths show poor correlations with H/3 line width. Table |4] 
lists the Spearman rank-order coefficients of these correla- 
tions. Since we found using FWTM and FWQM does not 
improve the correlations, we will focus on FWHM from now 
on. 

These results suggest that Civ and Cm] have different 
kinematics from Mg ll and the Balmer lines, and possi- 
bly originate from a different region than the low-ionization 
lines. However, since some dispersion in the Civ and 
Cm] FWHM is driven by the blueshift (e.g.. Fig. |4|, ac- 
counting for this dependence may reduce the difference 
in FWHM between Civ/Ciii] and the low-ionization lines. 
To test this, we plot the difference in FWHM Afwhm = 
log(FWHMH0/FWHMciii].civ), as a function of the blueshift 
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TABLE 4 

Spearman Test Results 





J-J 


^ran 


with FWHMh,3 






FWHMciv 


0.11 


0.39 


FWHMciii] 


0.14 


0.29 


FWHMMgii 


0.64 


2.8 X 10"" 


FWHMhc 


0.78 


1.8 X 10"'3 


FWHMciV.corr 


0.49 


6.7 X 10-5 


FWHMciIIJ.cor, 


0.36 


5.2 X 10-3 


with FWHMciv 






"'"' FWHMh;3 






logLi350 


0.02 


0.89 


log is 100 


-0.03 


0.83 


log(il350/i5100) 


0.21 


0.11 


EWCIV 


-0.07 


0.59 


K:iV-H/3 


0.59 


8.9 X lO-'' 


VciV-AlIII 


0.36 


5.2 X 10-3 


CIV AS 


-0.28 


0.028 



Note. — r, is the Spearman rank- 
order coefficient and Pran is the probabil- 
ity of being drawn from random distribu- 
tions. FWHMciv.corr and FWHMcmj.corr 
are the corrected FWHMs using the linear 
regression fits shown in Fig.|6] 

with respect to H/?for Cm] and Civ, in Fig.|6](top panels). The 
green lines are the linear regression results using the Baye sian 
method of Kelly (2 007) . We then use these best-fit linear re- 
lations to correct for the observed Cm] and Civ FWHMs: 

log FWHMciui/civ.con- = log FWHMciii]/civ + a + /3 AV , (2) 

where AV = (Voff_ciii]/civ ~ ^off.H/s) is the blueshift relative 
to H/3, and a and (3 are the best-fit coefficients of the lin- 
ear regression shown in green lines in the top panels of 
Fig. |6l [a,/3] = [0.062,-1.87 x lO'^] for Cm], and [a,/3] = 
[0.136,-1.57 X 10-4] for Civ. The "corrected" Cm] and Civ 
FWHMs are plotted against the broad H/3 FWHMs in the bot- 
tom panels of Fig. |6] This time significant correlations are 
detected for both Cm] and Civ, and Civ has the most sig- 
nificant improvement (see Table |4] for Spearman rank-order 
coefficients). Nevertheless, there is still substantial scatter 
(^ 0.15 dex for Cm] and ^ 0.12 dex for Civ) among these 
correlations. This "irreducible" scatter probably again reflects 
the different origins of Cm] and Civ from the low-ionization 
lines, which makes it difficult to bring their line widths into 
good agreement. We will return to this point in §15.11 

Although the blueshift relative to H/3 seems a viable proxy 
to correct the Cm] and Civ FWHM to better agree with the 
H/3 FWHM, it is of little practical value. One would like a 
proxy that can be determined from regions around the Cm] 
or Civ line alone. We have tried to correlate Afwhm with 
logLi35o, EWciv, asymmetry parameters of Civ (defined as 

AS = In I -^ I /111 ( T I ' where Ao is the peak flux wave- 

length, and A,ed and Abiue are the wavelengths at half peak flux 
from the model fits), and blueshift relative to Aim. We found 
that ApwHM is best correlated with asymmetry parameters of 
Civ and the blueshift relative to Aim at the Pran < 10"^ level, 
although still worse than the ones against the blueshift rela- 
tive to H/3 (where fian < lO"^). Using these weak correlations 
to correct for Civ FWHM does not seem to reduce the scatter 
between the Civ and H/3 FWHM much. 




-2000 2000 4000 

CIV-H/S blueshift (km s"^) 

4.2 




-2000 2000 4000 

Clll]-H/S blueshift (km s"^) 

Fig. 4. — FWHM as a function of the blueshift relative to the broad H/3 
centroid for CIV (upper) and CIII] (bottom). Both FWHMs seem to increase 
with the blueshift relative to H/3. Spearman correlation test results are r^ = 
0.44 (Pran = 5.1 X 10-^) for CIV and r, = 0.50 (Pran = 4.7 X 10-^) for CIII]. 

4.3. Comparing Single-Epoch Virial Mass Estimators 

The investigations so far in the previous two sections treated 
luminosity and line width independently. In principle, if there 
is covariance between line width and luminosity when com- 
paring the virial products based on different lines, the result- 
ing scatter in the residual virial products may be increased 
or reduced. Since we did not observe any strong dependence 
of line width on luminosity for any particular line, we expect 
such effects to be modest at most. 

As reasoned in the introduction, we adopt the H/3-I-L5100 
virial masses as the "truth" values, and minimize the 
differences using alternative line estimator with respect 
to H/3. There is more tha n one caUbration based 
on L5100 and FWHMh/? (eg.. iMcLure & DunlopI l2004t 
IVestergaard & Petersonll2006l lAssef et al.ll201 lb . and we use 
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Fig. 5. — Comparisons between different line widtlis and H/3 FWHM. Typical measurement uncertainties are indicated in each panel. The first row compares 
the FWHMs of Hq, Mg II, CIII] and CIV with H/3 FWHM. where the dotted lines show the unity relation. Only the FWHMs of Ha and Mg 11 show significant 
correlation with the FWHM of H/3. The second and third rows show similar comparisons, but with FWTMs and FWQMs for Ha, Mg II CIII] and CIV. The 
conclusion is the same. The bottom two panels show the co iTelations between FWTM/FWQM and FWHM for H/3. In each panel we show the best-fit linear 
regression results using the Bayesian method in IKellvl 120071 predicting Y at X): the best-fit slope, the uncertainty of the slope and the intrinsic random scatter 
about the regression. 

to be determined by linear regression analysis. 

In order to minimize the difference in virial masses com- 
pared to our fiducial masses, we allow the slopes on both lu- 
minosity and FWHM to vary. We u se the multi-d imensional 
Bayesian linear regression method in lKellyl ( l2007h to perform 
regression, treating our standard masses as the dependent vari- 
able Y, and (log L, log FWHM) as the 2-dimensional indepen- 
dent variable X. This approach takes into account measure- 
ment errors, and possible covariance between luminosity and 
FWHM, thus is arguably better than regressions on L versus 
L5100 and FWHM versus FWHMh^ separately. The regres- 



the IVestergaard & PetersonI (|2006[ VP06) calibration as our 
standard, which is compatible with our measurements of the 
H/3 FWHM, and provides simil ar estimates to those using the 
calibration in lAssef et al.l (12011, see below). The virial mass 
estimator based on a particular pair of line width and lumi- 
nosity is: 






1044 erg s- 



/FWHM\ 



(3) 
where L and FWHM are the continuum (or line) luminosity 
and width for the specific line, and coefficients a, b and c are 
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Fig. 6. — Top: FWHM ratio as a function of bluesliift, for CIII] (left) and CIV (right), compared with H/9. A significant congelation is detected for both CIII] 
and CIV. The green lines are the best linear regression fits using the Bayesian method in Kelly | 2()07 ). Bottom: "Corrected" FWHMs for CIII] (left) and CIV 
(right) using the best fits shown in the upper panels, compared with H/3 FWHM. A better one-to-one coiTelation is now seen between the CIIIJ/CIV FWHM and 
H/3 FWHM. although significant scatter still remains. 



sion results are listed in Table |5] 

In Fig.|7]we show the comparisons between different virial 
mass estimators and the VP06 H/3 estimator We show in the 
first two columns the comparisons for several calibrations in 
earlier work, and in the last two columns the comparisons for 
our new calibrations as summarized in Table |5] 

These earlier calibrations were calibrated using fainter sam- 
ples than probed here, but in general they provide mass es- 
timates that agree with the fiducial H/3-based virial masses 
(mean offset < 0.1 dex). The ex ceptions are the Mg ll-based 
calibrations in IShen et al.l ( 1201 lb . where the steeper slope of 
the 7?- L relation for Mg n than for Hf3 has led to increasingly 
larger discrepancies towards high luminosities^. For our new 



* We note that the Mg 11 calibrations in IShen et al.l COllj) were not based 
on linear regression fits against H/3 masses, and had a slope in the R — L 
relation fixed to be the one in McLure & Dunlop 1 2004). Using a steeper 
slope 0.62 in S-L3000 relation, the Mgll virial masses in lShen et alj J201 ih 



calibrations, the slope on FWHM is close to (albeit slightly 
smaller than) 2 for Ha, H/3 and Mg n, indicating that using 
FWHM in these calibrations improves the agreement with our 
standard mass estimator (VP06-H/?). However, for Civ, our 
linear regression result has a slope on FWHM that is much 
shallower This is because Civ FWHM is poorly correlated 
with H/3 FWHM for our sample, and the scatter between the 
two FWHMs rather than between the two continuum lumi- 
nosities is the dominant source of the diff erence in their virial 
masses (in contrast to lAssef et all 12011. see discussions in 
^5.11 ); therefore the regression prefers a smaller dependence 
on Civ FWHM to minimize the difference in the two mass 
estimates. In other words, the individual Civ FWHM adds 
little to improve the agreement with our standard mass esti- 

have negligible systematic offset relative to both H/?-based masses at z < 0.9 
and CIV-based masses at z > 1.5. 
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Fig. 7. — Comparisons between different virial mass estimators with tlie H/3 estimator in IVestergaard & PetersonI II2006I) . The left two c o lumns show 
the com parisons for existing calibrations based on either continuum luminosity or line luminosity, from l Assef et al.1 J201 IL Al 11 fShen et all j201 IL Sll), 
IGreene & Ho (2005, GH05), and Vestergaard & Peterson (2006, VP06). The right two columns show the results for our new calibrations based on our sample, 
where we allow the slope on FWHM to vary for altemative lines to minimize the residual in virial mass estimates when compared to H/3. The linear regression 
results are listed in Table [3] 

mates, but instead degrades the agreement. We found similar 
trends for Cm] (not shown) as for Civ. 

We test the dependence of the Civ virial mass residual on 
continuum luminosity, color and line shifts. The mass residual 
is correlated with log(Li35o/L5()(K)), but this is expected since 
both masses involve luminosity. Correcting for this color de- 
pendence only marginally improves the agreement between 
the two masses. On the other hand, the dependence on Civ- 
H/3 blueshift is strong enough such that incorporating this 
dependence can improve the agreement between Civ masses 
and the standard masses. These results again reflect the fact 
that the difference in FWHM is the dominant source in the 
virial mass difference. But this correction based on the Civ- 
H/3 blueshift is of little practical use since there is no need to 
correct Civ-based virial masses if we have H/3 coverage. Us- 
ing the Civ-Aliii blueshift or Civ asymmetry as a surrogate 
for the Civ-H/3 blueshift only leads to marginal improvement 
of the Civ-based masses, and thus is not of much practical 
value either 

Our new calibrations for Mg II yield consistent mass es- 
timates as those estimated from H/3 for luminous (Lsioo > 



10'*^'* ergs"') quasars. However, it would be useful to de- 
rive a Mgll calibration that is also applicable to lower lu- 
minosities. For this purpose we select ^ 900 z < 0.89 
quasars from the compilation in IShen et aTl (1201 Ih with 
good H/3 and Mgll measurements (AMvir < 0.1 dex) and 
LsiQo > 10^^ ergs"' (to reduce host contamination). We com- 
bine these quasars with the 60 high-luminosity quasars in 
our sample and perform the two-dimensional linear regres- 
sion on (logL30()o,logFWHMMgii) against the fiducial H/3- 
based masses. The best-fit coefficients are listed in Table 
H This Mgll calibr a tion i s close to the one presented in 
Vestergaard & Osmeij (l2009b . and works reasonably well for 



the entire luminosity regime lO''^ < L5100 < 10^' ergs 



5. DISCUSSION 

5.1. Comparison with Earlier Studies 

There have been many studies comparing different 
virial mass estimators (e.g., I Vestergaard & Peterson 2006| 
McGill et al. 2008; Dietrich & Hamann 2004; Dietric h et all 
2009- .Netzer et al...2007.; Shen et al...2008. .2011.; Wang et all 
|2009|) . These comparison studies used samples that have dif- 
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TABLE 5 








VirialMass Calibrations 






Type 


a 


h 


c 


offset 
(dex) 


(dex) 


Ref 


Fiducial mass 














FWHMhs,L5100 


0.91 


0.5 


2 


- 


- 


VP06 


Previous calibrations 














FWHMh<.,L5100 


0.774 


0.520 


2.06 


0.01 


0.14 


All 


FWHMH<.,Z-Ha 


1.239 


0.430 


2.10 


-0.08 


0.15 


Sll 


FWHMH,3,i5100 


0.895 


0.520 


2.00 


0.03 


0.003 


All 


FWHMH,3,iH/3 


1.676 


0.560 


2.00 


-0.04 


0.07 


GH05 


FWHMMgii,Z,3ooo 


0.860 


0.500 


2.00 


0.02 


0.16 


VO09* 


FWHMMgii,Z,3000 


0.740 


0.620 


2.00 


0.17 


0.17 


Sll* 


FWHMMgii,Z,Mgll 


1.933 


0.630 


2.00 


0.22 


0.16 


Sll* 


FWHMciv,Z.i350 


0.660 


0.530 


2.00 


0.10 


0.40 


VP06 


FWHMciv,iciv 


1.525 


0.457 


2.00 


0.09 


0.38 


Sll 


This work 














(L5ioo>10'*^'''*ergs-') 














FWHMhc«,L5100 


1.390 


0.555 


1.873 


0.01 


0.12 


_ 


FWHMhc<,Lhc 


2.216 


0.564 


1.821 


0.008 


0.15 


- 


FWHMh^,Lh,3 


1.963 


0.401 


1.959 


0.02 


0.04 


- 


FWHMMgii,Z.3000 


1.816 


0.584 


1.712 


0.02 


0.16 


- 


FWHMMgii,Z,Mgll 


3.979 


0.698 


1.382 


0.03 


0.15 


- 


FWHMciv,in50 


7.295 


0.471 


0.242 


0.03 


0.28 


- 


FWHMciv,i.civ 


7.535 


0.639 


0.319 


0.002 


0.26 
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Note. — Virial BH mass calibrations of Eq. (5] based on different line 
width and luminosity combinations for the 60 objects in our sample, calibrated 
against the VP06 FWHM-based H/3 virial masses. References of previou s cal - 
ibratio ns: GH05 JG reene & Ho 2005), VP06 ( Vestergaard & Peterson 2 00^ 
VO09 <Vestergaard & Osmen i2009i) . Sll (iShen et al.ii201 li) . All (Assef eTSl 
1201 11) . For each calibration we measure the average offset and scatter a rel- 
ative to the fiducial mass logMfu by fitting a Gaussian to the mass residual 
A = logM — logMflj. The scatter in the mass residual is dominated by intrin- 
sic scatter than by measurement eiTors. *In using the Mgll calibrations in 
[Shen et al. l20n|) and |%stergaard & Osmer (2009) we have added 0.125 dex 
to the measured logLsooo 2nd subtracted 0.067 dex from the measured logLMgii 
to account for the different fitting recipe used in this work. Note that th e absolute 
uncer tainties of these calibrations are on the level of > 0.3 dex (e.g.. lPetersoi]| 

Hon. 



ferent sizes, spectral quality, luminosity and redshift ranges, 
and focused on different lines. The current work is among 
the few studies that simultaneously investigate Civ through 
Ha in the same objects, and our sample is substantially larger 
and more homogen eous than used in simila r studies (e.g., 
[Dietrich & Hamannll2004tlDietrich et al.ll2009b . 

Our results agree with earlier work that t he line width of 
Mgii is well correlated with that of H/3 (e.g . , ISalviander et al.l 
l2007t lMcGill"etani2008l: IShen et al.ll2008f) . but now we have 
extended this conclusion to higher luminosity than can be 
probed with earlier samples at lower redshift. On the other 
hand, we confirm the poor correlation between Civ FWHM 
and H/? FWHM r eported earlier (e.g., iBaskin & Laorll2005l 
iNetzer et alJ 120071) . which was also inferred from the com- 
parison between Civ and Mg II using SDSS quasars (e.g., 
IShen et al. 2008). 

lAssef et al.. (.2011.) used a sample of high-redshift quasars 
with optical (covering Civ) and near-lR (covering the 
Balmer lines) spectroscopy to show that there is a corre- 
lation between the widths of Civ and the Balmer lines. 
They further concluded that the correlation persists, al- 
thou gh it becomes weaker, when other objects comp iled 
from lVestergaard & PetersonI (l2006b . iNetzer et al.l (l2007ll and 



[Dietrich et all (l2009h are included. 

The lAssef et al.l (1201 Ih sample is small (only 9 objects with 
all the measurements available for Spearman tests), and they 
do not probe a large dynamic range in Civ blueshift (see their 
fig. 13). The lack of large Civ blueshift objects in their sam- 
ple probably explains why they did not detect a significant 
correlation between the Civ virial mass residual and the Civ 
blueshift. 

To further investigate the disagreement on the correlation 
between Civ and H/3 FWHMs , we collec t ed lum inosity and 
FWH M measurements from i Assef et al.l (1201 ll 9 objects; 
Ain. I Vesterga^fe PetersonI (120061 21 o bjects; VP06) 
INetzer et al.l (l2007l 15 objects: N07), and iDietrich et al.l 
(200S| 9 objects; D09). N Ve use the P r escrip tion A measure- 
ments of Civ FWHM in lAssef etal ." (2011). Fig. [8] shows 
their distribution in the luminosity-FWHM space. The VP06 
sample probes a much lower luminosity regime than the other 
high-redshift samples. In the bottom panel of Fig. [8]we also 
show the ratio of FWHMciv/FWHMh^ against the contin- 
uum luminosity ratio L1350/L5100 (optical-UV color). The 
All sample spreads over a larger range in continuum color 
than the other samples, including se veral objects that ar e 
much redder than typical quasars (e.g.. [Richards et aTll2003l) . 
Only the D09 and All samples show a mild correlation be- 
tween FWHM ratio and optical-UV color, with Spearman 
rank-order coefficients of 0.65 (Pian = 0.02) and 0.67 (P,an = 
0.05), respectively. This correlation helps to reduce the virial 
mass differe nces between Civ and H/3 once this color effect 
is taken out (lAssef etal .11201 Ih . 

In Fig. [9] we show the comparison between Civ and H/3 
FWHMs for different samples. We run Spearman tests for 
the combined sample and for each subsample, and find that 
th e correlation betwe en the Civ and H/3 FWHMs reported 
in lAssef et aH (1201 ll) is essentially driven by objects in the 
VP06 sample, which probes a much fainter luminosity than 
other samples as indicated in Fig. [8] None of the other sam- 
ples show significant correlations between the two FWHMs, 
and this result does not change when we restrict ourselves 
to high-quality measurements. These high-redshift samples 
have a narrower dynamic range in line width than the VP06 
sample, and thus the intrinsic scatter between the Civ and H/3 
FWHMs can easily wash out any weak correlation. We there- 
fore reinforced our earlier conclusion in § |4.2| that. at least for 
the high-luminosity objects, the Civ FWHM is poorly corre- 
lated with the H/3 FWHM. 

5.2. Implications for High-Redshift Quasars 

The comparisons between different line estimators in pre- 
vious sections suggest that in the absence of Balmer lines, 
the Mgll estimator can be used as a substitute, which will 
yield consistent virial mass estimates to those based on the 
Balmer lines. On the other hand. Civ and Cm] can be used, 
although using the individually measured Civ/Ciii] FWHM 
does not seem to offer much advantage over simply using 
a constant value; of course, this conclusion is valid for the 
high-luminosity regime probed by this study. Civ is a compli- 
cated line, and may be more affected by a no n- virial compo- 
nent a s luminosity increases (see discussions in lRichards et al.l 
[201 ll) . These unusual properties of Civ suggest that it is 
likely the least reliable virial mass estimator at high-redshift, 
thus optical/near-lR coverage of Mg ll or the Bal mer lines is 
desire d for reliable virial mass estimates (e. g., Netze r et"an 
■2007;iTrakhtenbrot et al. 11201 ll:[Marziani & Sulentic..201 lb . 

One should also be aware that even for the most rehable 
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Fig. 8. — Comparisons between our sample and other samples in the lit- 
erature lAssefetal. (2011, 9 objects; All), Vestergaard & Peterson |2003 
21 objects; VP06), Netzer et al. (2007, 15 objects; N07), and Dietrich et 4\ 
{^09, 9 objects; D09)]. Upper. Distribution of H^ FWHM versus L5100. 
Middle: Distribution of CIV FWHM versus Lnso- Bottom: Distribution of 
the ratio FWHMciv/FWHMhs versus continuum color Li35o/i5ioo- Note 
that the VP06 sample has a much fainter luminosity than the other samples. 
This low-redshift sample also has smaller H/3 and CIV FWHMs than the 
other high-redshift samples. Only for the samples in Dietrich et al. ( 2003, 
9 object) and Assef et al. (2011, 9 objects) is there a significant correlation 
between FWHMciv/FWHMh;3 and Luso/Lsioo, which helps to reduce the 
virial mass difference between CIV and H/9 once this color-dependence is 
taken out. 
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Fig. 9.— Comparison between CIV FWHM and H/3 FWHM for different 
samples. The dotted line is the unity relation. Only for the low-redshift and 
low-luminosity VP06 sample is there a significant correlation between the 
two FWHMs. 

H/3-based virial mass estimates, there is still considerable 
scatter between v irial masses an d true masses (on the level 
of > 0.3 dex, e.g. lPetersonll201 ij) . Rare objects with unusual 
continuum and/or emis sion line properties (f or instance, dust 
reddened quasars, e.g., iRichards et al.ll2003h will lead to ad- 
ditional uncertainty in these virial mass estimates, along with 
measurement errors from poor spectral quality. Recognizing 
and accounting for the uncertainties in these virial mass es- 
timat es is crucial in essentially all BH mass related studies 
(e.g..'Shen et"ani2008l:lKe"llv et aUlIOOllMotlShen & KellM 
I20iai20ia) . 

6. SUMMARY 

In this paper we have empirically determined the relations 
between single-epoch virial mass estimators based on differ- 
ent lines for luminous (L5100 > 10'*^'* ergs"') quasars, using 
a sample of 60 intermediate-redshift quasars with complete 
coverage from Civ through Ha with good optical and near-IR 
spectroscopy. Our sample consists of typical quasars with no 
peculiarities in their continuum and emission line properties, 
has negligible contamination from host starlight, and is large 
enough to draw statistically significant conclusions. The main 
conclusions of this paper are the following: 

• The Mg n FWHM is well correlated with the FWHM 
of the Balmer lines up to high luminosities (L5100 > 
10'*^"' ergs"'), which justifies the usage of Mgn (in 
combination with L3000 or Livigii) to estimate virial BH 
masses for luminous quasars at high redshift. 

• The narrow-line contribution to the Civ line is gen- 
erally negligible for high-luminosity quasars; and the 
FWHMs of Civ and Cm] are well correlated, suggest- 
ing that both lines originate from similar regions. 

• The FWHM of Civ is poorly correlated with that of 
the Balmer lines, suggesting different BLRs for the Civ 
line and for the Balmer lines. Part of the discrepancy 
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between the Civ and H/3 FWHMs is correlated with the 
blue shift of Civ relative to H/3. 

• Using the FWHM of Civ increases the scatter between 
Civ and H/3 based virial masses, which is at least true 
for the high luminosity regime probed in this study. 
Cm] does not seem to be a superior substitute for Civ, 
both because of the non-correlation between Cm] and 
H^ FWHIVIs, and because Cm] is blended with Sim] 
and Aim. 

While in this work we focused on empirical relations, the 
correlations (and lack thereof) among these broad lines and 
continuum luminosities are ultimately determined by the ac- 
cretion disk and BLR physics. In future work, we will use the 
same sample to investigate in more detail the emission line 
and continuum properties of intermediate-redshift quasars, 
such as the Baldwin effect (Baldwin 1977), line centroid 
shifts, as well as average properties and correlations therein, 
in order to understand the underlying physics. In the mean 
time, we will expand our sample to include less luminous ob- 
jects and test these correlations over a larger dynamic range 
in quasar luminosity. 
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